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ABSTRACT 


Blood type testing was used to investigate the genetic 
makeup of Great Basin wild (feral) horse populations and 
answer the following questions: a) What is the extent of 
genetic diversity among wild horses? b) What is the relation- 
ship of wild horses to various breeds of domestic horses? c) 
What percent of foals within wild horse harem bands are sired 
by the dominant (alpha-male) stallion? Red blood cells and 
serum from peripheral blood samples of nearly 1000 wild 
horses from five areas were tested by serological and 
electrophoretic techniques for genetic markers at 19 polymor- 
phic” 1ocis The average number of alleles for 13 domestic 
breeds was 73.3+8.0 (range 58 to 86). The average for wild 
horse populations was 69.7+9.0 (range 57 to 82). The average 
frequency of heterozygotes per locus (average heterozygosity? 
averaged 0.353+0.011 (range 0.295 to 0.443) for domestic 
breeds. Average heterozygosity for wild horse populations 
averaged 0.402+0.009 (range 0.368 to 0.442). Dendrograms 
constructed using pairwise comparisons of Nei’s distance 
measurements (D) for the domestic breeds and the wild horse 
populations substantiate anecdotal accounts of the origins of 
Great Basin horses from draft horses, saddle horses of 
American breed origin and Spanish Barbs. Analysis of genetic 
markers showed that nearly one-third of foals were not sired 
by the harem stallion of the band into which the foal was 
born. Harem stallions averaged 12.2+5.1 years (range 7-25) 
and dams averaged 7.1+4.0 years (range 2-20), about a genera- 
tion younger than the stallions. 


I. INTRODUCTION 


Following intense public pressure to provide Federal 
protection for wild (feral) horses and burros from inhumane 
Slaughter and possible extinction, in 1971 the Wild and 
Free-Roaming Horse and Burro Act was passed. The bill 
provides for the protection, management, and control of wild 
horses and burros that roam public lands administered by the 
Department of Interior through the Bureau of Land Management 
(BLM) and by the Department of Agriculture through the Forest 
Service. — 


With no natural predators and under protection of the 
Federal law, the numbers of wild equines dramatically in- 
creased after passage of the 1971 act, requiring a plan for 
management and removal of excess animals on public lands. 

The Public Rangelands Improvement Act of 1978 (PRIA) called 
for a scientific committee to outline a research study that 
would further knowledge "...of wild horse and burro popula- 
tion dynamics..." to assist the Secretary of the Interior 
".se2e1M making his determination as to what constitutes excess 
animals." 


At the request of BLM, a committee of the National 
Research Council (NRC) issued a report outlining such re- 
search projects. Blood type testing was proposed to inves- 
tigate the genetic makeup of wild horse populations 
(Committee on Wild and Free-Roaming Horses and Burros, 1980). 
Such testing could answer the following questions and thereby 
provide a basis for management decisions: 

a. What is the extent of genetic diversity among wild 
horses? 

Dis What is the relationship of wild horses to various 
breeds of domestic horses? 

c. What percent of foals within wild horse harem bands 
are sired by the dominant (alpha-male) 
stallion? 


This report outlines our findings concerning WILD HORSE 
PARENTAGE AND POPULATION GENETICS conducted under the 
auspices of the United States Department of the Interior, 
Bureau of Land Management Solicitation AABS2-RPS-27. Genetic 
data from nearly 1000 wild horses sampled at seven trap sites 
in the Great Basin area of Oregon and Nevada were obtained by 
1) blood typing peripheral blood samples for markers at 19 
loci and 2) observation of coat color and pattern variants at 
seven loci. 


IIT. MATERIALS AND METHODS 
A. Study sites of wild horse populations 


With the advice of BLM horse specialists, five sites 
were chosen in the Great Basin area between the Sierra Nevada 
and Rocky Mountains of the western United States. Two sites 
had two trap locations each, providing paired subpopulations; 
the other three sites had one trap site each. The sites were 
located so far apart that the possibility of migration be- 
tween sites could be ignored. 


Sites were chosen to meet requirements of both the 
Parentage and Population Genetics project conducted by the 
University of California, Davis (UCD) and the Fertility 
Control projects of the University of Minnesota/Marshfield 
Medical Foundation (MM). Fertility control projects to as- 
sess methods of limiting population growth included stallion 
Sites to study the effectiveness of vasectomy for dominant 
males and mare sites to study the feasibility of hormone 
implants. At least 30 bands were needed for stallion project 
areas, and at least 100 breeding age females for the mare 
areas. Adequate sample numbers for genetic studies could be 
obtained within the requirements of the Fertility Control 
project. Since the Fertility Control studies would be con- 
tinuing for at least two years, sites were selected which 
would not be subject to BLM reduction roundups during that 
time. For paternity data, sites should not have been dis- 
turbed by roundup or animal removal subsequent to conception 
of the foals (gestation length = 11 months). All sites sup- 
ported populations of over 100 animals but 1n no case were 
all the horses captured. 


1. Flanigan (F) was located in northwestern Nevada, northwest 
of Pyramid Lake, Washoe County, north of Reno. In 
September 1985 BLM had removed over 300 horses in 
this area. In December 1985, we obtained blood 
samples from 175 horses in 3SO bands. Band composi- 
tion may have been affected by the recent BLM 
gather. 


2. Wassuk (W) was located in the Wassuk mountains of western 
Nevada, west of Walker Lake, Mineral/Lyon County, 
south of Yerington. Horses in this area had been 
the subject of a thesis by Pellegrini (1971) 
describing territoriality and movement patterns of 
wild horses. Reduction roundups had never been 
conducted by BLM at this site. The roundup and 
blood sampling of 119 horses in 21 bands was con- 
ducted in January 1986. 


3S. Beaty Butte (B) was located in extreme southeastern 
Oregon, Harney County. Horses in this area had 
been the subject of previous studies by Eberhardt, 
Majorowicz & Wilcox (1982) to determine reproduc- 
tive rates in feral horses. A BLM reduction roun- 
dup had been conducted in November 1984. Foals in 
Our gather of February 1986 had been born in spring 
1985 (conceived in spring 1984), and band structure 
relative to foal paternity considerations may have 
been affected. From this area 112 horses in 17 
‘bands were blood typed. 


4. Stone Cabin 1 (SC1) was located in west central Nevada, 
Nye county, east of Tonopah, north of highway 6. 
Horses in this vicinity had been described by Green 
& Green (1977) in a study of demography, habitat 
utilization, social organization and behavior. No 
reduction roundups which might influence the pater- 
nity studies had been conducted at this site. Our 
gQather of 239 horses in 30 bands was in September 
L7oo. 


S.- Stone Cabin 2 (SC2) was located about 10 miles from SCl, 
south of highway 6. Our gather in September 1986 
consisted of 127 horses which ran together as one 
band so this site generated no paternity data. 
Blood types could be used for gene frequency data. 


6. Clan Alpine 1 (CA1) was located in west central Nevada, 
Churchill County, between Fallon and Austin in the 
Clan Alpine mountains. No reduction roundups which 
might influence the paternity studies had been con- 
ducted at this site. Altogether 104 animals in 17 
bands were gathered in October 1986. 


7. Clan Alpine 2 (CA2) was located in the same area as CAI, 
about 10 miles away. Altogether 99 horses in 12 
bands were gathered in October 1986. 


Horses were gathered by helicopter using BLM-approved 
contractors, field methods and equipment. Two contractors 
were involved, one for the first five trap sites and a second 
for the two CA sites. The expertise and skill of the 
helicopter pilots was relied on to identify and gather intact 
bands on a band-by-band basis. An animal behaviorist was 
present at all but the Wassuk site to help identify dominant 
stallions and mare-foal pairs. Confined in restraining 
chutes, horses were were aged by evaluating tooth replacement 
and wear and bled by jugular venipuncture. Sex and color 
were recorded. Appropriately selected horses from all but 
the CA area were chosen for fertility studies and tagged with 
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identifying materials including numbered neckbands and radio 
transmitters before being released. Horses in the CA gather 
were not part of the Fertility Control study, and were 
removed from the area after capture. 


C. Domestic horses 

Blood samples were sent to the laboratory in conjunction 
with breed registry regulations requiring blood typing or re- 
search programs studying genetic diversity within and among 
domestic horse breeds. Thirteen breeds, their abbreviations 
and number of samples included in gene frequency determina- 
tions were: Arabian (AR) (5448), Thoroughbred (TB) (14517), 
Quarter Horse (QH) (1195), Standardbred (ST) (2959), Morgan 
(MH) (1006), American Saddlebred (SB) (803), Tennessee 
Walking Horse (TWH) (121), Belgian (BE) (771), Shire (SH) 
(138), Argentine Criollo (ARC) (121), Chilean Criollo (CHC) 
(194), Mangalarga (M) (120) and Mangalarga Marchador (MM) 
CSF 


Breeds were selected which have been historically and 
anecdotally associated with Great Basin wild horses. Each of 
these breeds is presently maintained with a studbook record. 
ARs and TBs are considered foundation stock for modern breeds 
of horses and were often used by the Army as a source of 
Remount stallions. Five of the breeds (QH, ST, MH, SB, TWH) 
have been developed in the US primarily as saddle horses from 
stock brought in by early settlers. Although the origins of 
these five breeds are similar, through selection for specific 
uses (e@.g-, racing, showing, cattle work) each breed now has 
distinctive characteristics. Draft breeds are represented by 
BE and SH. Horses whose origin is thought to be nearly ex- 
clusively from the Iberian peninsula during the 16th century, 
at about the time the Spanish settlements and missions of the 
Southwest were bringing horses to the US, are represented by 
four Andalusian/Barb South American breeds: Criollo horses 
from Argentina (ARC) and Chile (CHC) and two breeds from 
Brazil (M and MM). 


D. Blood collection 

Blood was collected by venipuncture into evacuated blood 
collecting tubes of approximately 10cc draw. One tube con- 
tained the anticoagulant acid-citrate-dextrose (ACD) from 
which red blood cells (RBCs) were obtained for serological 
and electrophoretic studies of blood type markers. The other 
tube contained no additives and was used aS a Serum source 
for additional electrophoretically determined markers. Blood 
samples were maintained on ice or under refrigeration until 
tests were completed. 


An extensive battery of tests for detecting genetic dif- 
ferences has been developed and applied to the pedigree 
verification of purebred horses (e.q., Braend, 1973; 
SemOi@m welt 7 eae Ott. Lo 7es Bowling & Clark, .1989). All 
loci selected for routine testing are polymorphic. The tests 
are relatively easy to perform, accurate and highly effective 
in recognizing incorrect parentage. It was assumed that 
these tests would be effective for answering questions of 
wild horse paternity. The large database available for 
Gomestic horses could be utilized for studying genetic 
relatedness of wild horses to domestic breeds. 


Nomenclature was in accord with internationally agreed 
usage, except for variants at some loci for which an interna- 
tional nomenclature has not been adopted. 


1. Serological testing for blood groups 


Standard immunological procedures involving hemag- 
Qlutination and complement mediated hemolysis (Stormont & 
Suzuki, 1964; Stormont, Suzuki & Rhode, 1964) were used to 
detect red cell alloantigens at seven internationally recog- 
nized blood group loci (Table I). 


a a a a a a a a a ae ae eae a 
cm a a a a ae ae ae ae 


Table I. Serologically detected RBC alloantigenic factors 


Genetic Allo-antigenic Determinants Number 
Systems (Reagents) of Alleles 
A Aa,Ab,Ac,Ad,Ae,AF,AQG iz 
Cc Ca Z 
D Da, Db, De, Dd, De, DE, Dag, Dh, Di, Dk, 24 


D1, Dm, Dn, Do, Dp, NF*12,NF19,NF23, 
NF44,NF69,NF71,NF98, NFO8 
K Ka Ps 
e Pa,Pb et 
Q Qa, Q@b,Qc = 
U Ua 2 
Total: 530 
XNF: laboratory designation for new factors, not yet 
internationally recognized. 
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An RBC sample from each horse was mixed with an aliquot 
of each specific reagent. A positive reaction was recognized 
either as agglutination (clumping of the cells) or hemolysis 
(rupture of the cells with the addition of rabbit serum as a 
complement source). Positive reactions detected determinants 
comprising phenogroups (alleles) whose inheritance has been 
well documented in domestic horses. A genetic type was as- 
Signed for each animal based on the reaction pattern of its 
RBCs to the battery of reagents. 


2. Electrophoretic testing of serum and red cell enzymes and 
other proteins 


Electrophoresis and isoelectric focusing were used to 
analyze biochemical differences in blood constituents 
(Braend, 1974; Sandberg, 1974; Juneja, Gahne & Sandberg, 


1978; Braend & Johansen, 1983). RBC lysates or serum samples 
were inserted in a starch or polyacrylamide gel matrix, then 
an electrical current was applied. When standards included 


in each gel migrated an established distance, the current was 
shut off and the gel stained with protein dyes or biochemical 
reagents. Genetic variants in each sample were named in com- 
parison with known reference samples within the same gel. 
Twelve systems of proteins were analyzed by means of this 
technique (Table II). 


nn iii 


Table II. Electrophoretically detected genetic variants 


Systems Number of Alleles 
Albumin (Al) ae 
Transferrin (Tf) 14 
Xk 4 
Esterase (Es) ry. 
Protease Inhibitor (Pi) 24 
Gc 2 
Carbonic Anhydrase (CA) 6 
Catalase (Cat) ve 
Phosphoglucomutase (PGM) > 
Phosphohexose Isomerase (PHI) 4 
Phosphogluconate Dehydrogenase (PGD) 3 
Hemoglobin-alpha (Hb) a 

Total: 76 


me ee ee me we me a a a wa ae a ae ae 
eee ee me ee er ee a a aa 
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F. Gene frequencies, heterozygosity and genetic distances 


Blood group frequencies were calculated as follows: for 
C, K and U by the square-root method assuming Hardy-Weinberg 
(H-W) proportions; for P, with Bernstein’s equation developed 
for human ABO; and for A and Q, with the square-root method 
to obtain frequencies of the null alleles, followed by cal- 
culation and adjustment of the other allele frequencies to 
obtain the best fit to the data under the assumption of H-W 
equilibrium; for D, by direct counting assuming no ambiguous 
phenotypes. Allelic frequencies for protein polymorphisms 
were determined by direct counting from phenotypes. 


Computer programs of K. Ritland (Department of Botany, 
University of Toronto, Toronto, Ontario, Canada) were used to 
calculate average heterozygosity, Nei’s measures of standard 
genetic distance (D), total gene diversity or heterozygosity 
(Hy), diversity within subpopulations (breeds or trap sites) 
(Hs), diversity among subpopulations (Dco7), diversity of sub- 
populations relative to the total population (GsrT) and to 
construct a dendrogram from Nei’s genetic distances using the 
unweighted group pair method. The number of effective al- 
leles (Hedrick, 1983) was calculated from average 
heterozygosity. 


Gene frequencies for nine coat color genes, each with 
two alleles, were calculated for wild horses, assuming H-W 
proportions, by sequentially sub-dividing each population 
(Trommershausen-Smith, 1979). First, pink-skinned white 
horses were counted as heterozygous for W and then removed 
from further consideration because epistasis makes it vir- 
tually impossible to ascertain their alleles at other coat 
color loci. Second, the frequency of the gene controlling 
gray color (G) was one minus the square root of the propor- 
tion of non-grays in the population. The grays were not fur- 
ther considered due to the epistatic interaction of G with 
the other coat color genes. Among the remaining horses, roan 
(R), overo (including sabino) (0), dun (D) and tobiano (T) 
were infrequent in all populations and were directly counted 
as due to the action of the respective dominant genes in 
heterozygous condition. The gene frequency for C at the al- 
bino locus was one minus the proportion in the non-gray 
population of palominos and buckskins. At the extension 
locus, E was one minus the square root of the proportion of 
red horses (chestnuts, sorrels, red duns, palominos) in the 
non-gray population. At the agouti locus, A was one minus 
the square root of the proportion of blacks and grullas to 
the aggregate of blacks, bays, browns, buckskins, and 
grullas. 
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G. Parentage determination 

Named blood type factors behave as if controlled by 
codominant genes, and epistatic interactions between loci are 
not known to occur. Null alleles for which no blood group 
antisera are available or no protein product is discernable 
upon electrophoresis behave as recessives to named factors. 
Two rules of Mendelian inheritance are invoked in the assess- 
ment of parentage. According to Mendel’s first law, if a 
factor is present in an offspring, at least one parent must 
also have that factor (principle of dominance). According to 
Mendel’s second law, one member of each allelic pair must be 
transmitted to each offspring (principle of segregation). 


Applying the principle of segregation, a foal was ex- 
cluded as an offspring of a particular sire if it failed to 
share with the stallion an allelic specificity in one or more 
systems. If no incompatibility could be determined, the foal 
was scored as qualifying as an offspring of the stallion. 

The effectiveness of blood typing to detect incorrectly as- 
Signed paternity was at least 98% (Bowling, 1985) if a dam 
was also identified for a foal (allowing application of the 
Principle of dominance to the parentage determination.) 
Mare-foal pairs were difficult to ascertain at F, W and B due 
to the near-weaning age of the foals, but pairings were re- 
corded for the SC1 and CA sites. Usually blood typing could 
be used to assign a dam to a foal within a band, even without 
benefit of observation records. For F, W and B, ina total 
of three cases a mare within a foal’s band could be found 
that qualified as its dam, but the foal was excluded as a 
product of the mating of this mare and the dominant stallion. 
In such cases the likelihood of the foal being associated 
with its true dam seemed greater than with its sire, so the 
dam relationship was counted as correct. The sire was ex- 
Cluded even though he could not have been excluded without 
considering a dam. 


For this study, a qualification was interpreted as proof 
of parentage, although it would not constitute legal proof. 


III. RESULTS AND DISCUSSION 


Genetic variation provides the flexibility for an in- 
dividual to adapt to diverse environments and for populations 
to have evolutionary potential to survive in the face of en- 
vironmental changes (Dobzhansky, et al., 1977). The process 
of mutation is the source of variation, but it is very slow 
and most new mutations are lost by random drift within a few 
generations of their appearance. Genetic variation within a 
population is thus an extremely important reservoir which 
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must be protected to ensure the welfare and evolutionary 
potential of the species. 


1. Gene frequencies and H-W equilibrium proportions 


All blood samples were received at the laboratory 
in excellent condition for testing. Gene frequencies for the 
seven trap sites of Great Basin wild horses are presented in 
Appendices 1A-G. Among the 19 loci in the wild horses, three 
sites were fixed at the K locus and one at Gc. To calculate 
gene frequencies for blood group loci which have null al- 
leles, each trap site was assumed to consist of animals from 
a random-mating population for which observed and expected 
phenotypes would not show significant deviation from H-W 
proportions. This assumption was substantiated using protein 
loci for which genotypes could be directly inferred from 
phenotypes. Chi-square tests of H-W proportions showed no 
Significant deviation for any trap site at the simple Al 
locus in which two alleles each usually occur in moderate 
frequency. A few minor deviations seen at other simple loci 
(Xk, Ge, PGD, Hb, Cat) were not significant when classes were 
combined to eliminate classes with small numbers. Thus we 
verified the assumption that horses sampled at each trap site 
could be considered as members of an interbreeding 
population. Furthermore, the H-W proportions indicated that 
assumptions of diagnosis and inheritance of variants based on 
studies in domestic breeds were correct for wild horses. 


The five populations were geographically isolated from 
each other. However, horses from paired trap sites within 
each of two populations (SC and CA) could have been members 
of a single random-mating population. Combining the two 
sites of CA data as though from one population provided no 
Significant deviation from H-W proportions. Significant 
deviations from H-W proportions were observed if the two SC 
trap sites were combined into one data set. Both pairs of 
trap sites were only separated by 10-20 miles, but the 
presence of a highway between the two SC sites may have been 
the effective barrier which prevented the two sites from be- 
having as one population. 


2. Total alleles, effective alleles, average heterozygosity 


Comparison of allelic variation in horse populations is 
shown in Table III. Using a population size base of 100 (200 
alleles/locus), both uncommon (gene frequency between 0.05 
and 0.01) and rare (frequency between 0.01 and 0.005) alleles 
are included as well as the more common polymorphic alleles 
with frequencies greater than 0.05. Also given in Table III, 
calculated from gene frequencies, are values for the effec-— 
tive number of alleles and average heterozygosity. 
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Table III. Allelic variation in horses at 19 blood type loci 


Trap No. of Alleles Nov’ OF Average 
Seve. (8-2 ————F esse => S-=- Effective Hetero- 
Breed big Crp Protein’ Total Alleles zyQgusity (SE) 
F 3O 41 wy 55.3 OF 5/E"tO., 0585 
W 24 SS D7 Servs GO. 580° (OV062) 
B 2a 36 61 id 0.442 (0.045) 
Sct 30 48 78 46.3 OF4823">'(O. 062) 
SEZ 28 38 66 43.0 Oo. 4769407059) 
CA1 S1 4 B2 42.5 0.404 (0.058) 
CA2 28 45 vs 40.5 O07 368 "'(0,'061 } 
AR 29 S/ 68 oe | 0.346 (0.062) 
TB 24 34 58 SS. 7 OFZ7TO= {OFO6O) 
QH 36 30 B& 45.0 0.403 (0.062) 
Si 29 40 69 OTE tT 0.413 (0.048) 
MH 34 ow 85 42.9 OF410>(OF056) 
SB 34 47 81 40.7 OF 2586.0. 059) 
TWH 2? 42 wi So. 67550 €0. 056) 
BE 31 42 7S 46.8 0.443 (0.059) 
SH PAZ] 43 Ue SB. O7361"7°0.058) 
ARC 27 46 7 Ae 41.1 OF47O* (07061) 
GHC 30 ag 74 41.8 OF 4267 (0.052) 
M 26 Se 64 26.3 Ovo0o) XO. 068) 
MM 3O 48 78 45.4 0.412 (0.066) 


The total mumber of alleles detected by tests for 19 
blood type loci was 126, but no breed or population contained 
them all. The average number of alleles for domestic breeds 
was 73.3+8.0 (range 58 to 86). The average for wild horse 
populations was 69.7+9.0 (range S57 to 82). 


The number of effective alleles is an inverse function 
of the theoretical homozygosity of the population. This 
measure provides a comparison of populations for number of 
alleles, while minimizing the effect of infrequent, uncommon 
or rare alleles. The number of effective alleles for domes- 
tic breeds averaged 40.3+4.0 (range 33.7 to 46.8) and for 
wild horses averaged 41.3+2.8 (range 38.8 to 46.5). 


The average frequency of heterozygotes per locus 
(average heterozygosity) considering the 19 polymorphic loci 
averaged 0.353+0.011 (range 0.295 to 0.443) for 13 domestic 
breeds. Average heterozygosity for wild horse populations 
averaged 0.402+0.009 (range 0.368 to 0.442). All the wild 
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horse populations had higher heterozygosity values than the 
average for domestic breeds maintained within studbooks, al- 
though the highest value was not greater than found in the 
domestic breeds. 


Many domestic horse breeds have breeding practices which 
would tend the restrict gene pools and decrease 
heterozygosity. Some breeds (e. g., AR, TB, MH) are main- 
tained with a closed studbook which does not allow the 
registration of animals whose dams and sires are not in that 
studbook or in related studbooks from other countries. New 
Variants within closed studbooks would be due only to un- 
detected introgression or to mutational events. In some 
breeds linebreeding or inbreeding which would tend to 
decrease heterozygosity are used as a tool to intensify the 
traits of highly valued horses. Despite these practices con- 
Siderable amounts of genetic variation are found within 
domestic breeds and average heterozygosity values for the 
polymorphic loci of this study do not vary dramatically. 


Considerable genetic variation was found in the Great 
Basin wild horse populations just as in breeds of domestic 
horses. This variation suggests that the origin of each 
population probably was not from a severely restricted foun- 
der population or that the population did not pass through a 
bottleneck. Although examples of inbreeding have been 
reported in wild horses (Kirkpatrick & Turner, 1986), Duncan 
et al. (1984) showed effective avoidance of inbreeding for a 
small group of Camarque wild horses. Neither the allelic 
data nor the heterozygosity data provided substantial 
evidence in these Great Basin wild horse populations for 
restricted gene pools which could be a consequence of 
inbreeding. 


Heterozygosity has also been found to be characteristic 
of another equine species, the true wild horse. Przewalski’s 
horses (Equus przewalskii) exist only in captive populations. 
Based on 13 founders about 11 generations ago, the species is 
highly inbred, but shows a surprisingly high level of 
heterozygosity for a restricted number of blood type variants 


(Bowling & Ryder, 1987). 


Additional population measures (Nei, 1973) are shown in 
Table IV. Hy is total heterozygosity, or the probability 
that two gametes chosen at random from the total population 
will carry different alleles. As with the related measure of 
average heterozygosity in Table III, wild and domestic horses 
do not differ greatly. Hs 15 subpopulation heterozygosity, 
which is nearly as large as Hy for any Table IV populational 
cluster. Interpopulational diversity (Ds7) is low, both in 
wild (0.027) and domestic (0.060) horses. The diversity of 
subpopulations relative to the total population (Gs7) was 
smaller for wild horses (0.063) than for domestic breeds 
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fO.2%36). Thus the differences between populations of Great 
Basin wild horses are less than the differences between 
breeds of domestic horses. 
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Table IV. Nei’s population measures for horses 


Population Cluster 


7 Wild + 13 Domestic 7 Wild 13 Domestic 
Hy 0.441 0.429 0.443 
Hs 0.3589 0.402 0.383 
Dst 0.052 0027 0.060 
Gst Ost 17 0.063 0.136 


3. Unique variants 


Except at the Pi locus, no unique variants were observed 
in wild horses which had not been found in domestic breeds. 
The unusual Pi types are recorded as "oth" in the gene fre- 
quency listings (Appendices 1A-G). The highly polymorphic Pi 
locus is thought to be a complex locus with more than one 
protein product (Juneja, Gahne & Sandberg, 1979; Pollit & 
Bell, 1983; Patterson & Bell, 1987). Genes within the com- 
plex may recombine to produce new protein combinations and 
thus show evidence of new variants more often than a trait 
produced by the product of a single locus. 


A striking association was seen between wild horses of 
SC1, CA1 and CA2 and the Iberian peninsula-derived breeds of 
M and MM for a fast-migrating transferrin variant designated 
Tf-A. In the Brazilian breeds, the frequencies were 0.008 
(M) and 0.165 (MM), and in the wild horses 0.04 (SC1), 0.014 
(CA1) and 0.126 (CAZ2). This variant has not been reported in 
scientific literature and outside the sites/breeds listed 
above has only been recognized in this laboratory on six 
previous occasions among about 150,000 samples. The sharing 
of this very rare variant may be due to a relatedness in 
Origin between those wild horse groups and horses from the 
Iberian peninsula. Alternatively, the occurrence of the 
variant in the two groups may be due to coincidental inde- 
pendent mutation for the same variant or for two variants 
which appear to be the same using standardized conditions of 
electrophoresis. 
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4. Coat colors 


Gene frequencies for coat colors are presented in Table 
V0 Due to the unavailability of coat color gene frequencies 
for domestic breeds and the great potential in wild horses 
for misassignment of color due to conditions of foal coat and 
long winter hair, this information was not used in calcula- 
tions of genetic distance. 


Trap Boat Color Loci. 

Site@l/  -r rr 
Breed W G R T D 0 C s A 

F On lO 2 me OLY mer O Om OO. 0 O.O00) O00. O.49 0.61 
W 0.00 rare* 0.06 0.00 0.00 0.00 0.99 0.54 0.80 
B reme  O.O4 weare.rare 0.02 0.00 0.99 0.68 0.64 
Sst oan SO eee Oe eee Orono.  O..OO. O..965 . OO, 357. 0.6) 
Siz Oey Cm me OO. Oe OOM OOO OOO 1,00, 0. 02 0.47 
CA1 OS Ore On Ot O.OO O02 rare .O,97 .0.33 0.97 
CA2 eee Wee Oa nO. OO. OF OO OSS Onze 60371 


Gene frequencies are given for the dominant allele at each 
locus. The frequency of alternative allele is one minus 
the given number. 

2Rare: frequency equal to or less than 0O.OOS. 
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Coat color variation was prominant in wild populations 
only at E and A, which is typical for domestic horses, as 


well. Several breeds express variants at G, R, D and C loci, 
but variation at W, O, T is generally confined to breeds 
which selectively breed for coat color variants. Dominant 


alleles at some of the loci (W, O, R) may be lethal in 
homozygous condition (Hutt, 19693; Hintz & Van Vleck, 1979; 
Vonderfecht, Bowling & Cohen, 1983) but are so infrequent in 
the wild populations that their effect on compromising 
population growth is probably negligible. 


Nei’s distances (D) between each wild horse population 
and the 13 domestic breeds are presented in Appendices 2A-G 
as bar charts for each trap site. Dendrograms constructed 
from genetic distance comparisons for various combinations of 
the 20 breeds/trap sites by the unweighted group pair method 
are shown in Appendices SA-E. 
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Nei’s distance is a measure of genetic relatedness be- 
tween two populations which may reflect the historical and 
phylogenetic record. Gene frequencies and thus D can be 
Grastically affected by mutation, natural selection, artifi- 
cial selection or geographical isolation and lead to spurious 
associations if not interpretated with caution. A sequential 
plotting of dendrograms for several arrays of population sub- 
sets is useful for discerning those associations between 
breeds which remain consistent, from others which may be ar- 
tificial and spurious. 


Dendrograms obtained from the blood typing data for nine 
US domestic breeds are consistent with breed origins. The 
first dendrogram (Appendix SA) shows three major clusters, 
one for the five American breeds (QH, ST, SB, MH and TWH), 
one for the draft horse pair (SH, BE) and one for the TB-AR 
pair. 


A second dendrogram (Appendix 3B) constructed of 13 
domestic breeds shows the South American breeds clustered 
with the US breeds, in opposition to both the draft breeds 
and the AR-TB pair. The clustering of TWH with the three 
South American breeds may be an artificial association due to 
selection to produce horses that consistently perform special 
Qaits. No direct connection exists between TWH and the 
breeds with which it is clustered; on the other hand the as- 
sociation with MM, which also is a gaited breed, may reflect 
a common connection to Iberian horses which performed at the 
amble or pace. 


The third and fourth dendrograms (Appendix 3C,3D) were 
constructed of the nine US breeds plus six Great Basin wild 
horse trap sites. Only one CA site was included in each 
dendrogram, because using both caused them to cluster 
together rather than pointing to their relationships with 
with other clusters. Four trap sites (SCZ, CA1,SC1 and F) 
formed a cluster which, with the addition of BE, appeared to 
be a split from four of the American breeds (QH, MH, SB and 
ST) plus CAZ. Two other wild horse trap sites (W and B) ap- 
peared to be distinctive from the other sites. 


A #ifth dendrogram (Appendix SE) consists of all the 
populations and subpopulations of this study. The South 
American breeds of Iberian origin primarily appear as add-ons 
to existing clusters and do not seem to have a primary role 
in the origin of either the wild horses or the other domestic 
breeds, although they are included within clusters for all 
the groups except for SH, TB and AR. 


With respect to wild horses these dendrograms are basi- 
cally consistent with Amaral’s (1977) hypothesis of the 
Origin of Great Basin wild horses from escaped or released 
draft horses, cavalry mounts, and saddle animals as well as 
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mustangs derived from Spanish Barb stock. The cavalry 
horses, often of AR or TB origin, seem to have had a minor 
role especially compared to American breeds of saddle horses. 


1. Social organization of wild horse populations 


Numerous studies (summarized by Kirkpatrick & Turner, 
1986) have described the family units (bands) which organize 
wild horse populations. A harem band consists of a single 
stallion, mares (generally 1-5) and their offspring. 
Multi-male harem bands have been reported in which usually 
only one stallion demonstrates dominance (alpha-male). Harem 
bands are generally described as having a stable composition, 
with changes in band composition usually the result of youn- 
ger animals leaving the herd, or death of the dominant 
stallion. All-male bachelor bands are another type of group 
found in wild horse populations, but generally these bands 
are less cohesive than harem bands. 


Table VI presents the age, sex and band data for the 
seven wild horse trap sites. The SC2 site was distinctive 
from the other sites in its loss of band integrity in 
response to the helicopter threat. 
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Table VI. Bands, sexes and ages for wild horse populations 


Trap Sites 


F W B SC1 SEZ CA1 CA2 
N of horses 1405 LS, Na ar PERS 12Z/ 104 99 
Bands 30 Zi i7 50 1 le] ig2 


Bancdweize canoe 2-15. S-13 1-10. $=14 NA 2-14 2-14 


N of foals 23 daz ta! 45 26 24 Fal 
Female 16 8B 20 18 1a NV? aul 
Male if, 4 } ie) Zi 10 ars 10 

N of yearlings ze 7p Ze Za het 6 6 
Female iS S ro) =) 4 1 4 
Male 8B me 16 Paik ys 5) ¥e 

N aged 2-5 64 oF ZS 88 49 ay Sil 
Female 34 20 ee, 49 Pa) 20 ie 
Male 3O yee? 4 oo 24 19 14 

N aged 6+ 67 61 2 dal 4} os 41 
Female Spe) 28 16 41 18 Ze 18 
Male Sz KG 16 36 Zs aS 23 
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2. Results of paternity determinations 


ee ) Field observations suggest most foals in wild horse 
populations are sired by the dominant stallion of the harem 
band into which the foal was born. Although occasional 
sightings suggest that subdominant stallions may breed mares, 
the extent to which foals are sired by stallions other than 
the harem male is not known. 


For the assessment of paternity our data included 69 in- 
tact harem bands with a total of 121 foals (Table VII). 
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Table VII. Foal paternity in harem bands 


Trap Site 


F 4 B SC1 CAL CA2 Total (%) Total® (4) 


HAREM BANDS 12 10 12 20 8 7 69 45 
BANDS, 1 stallion 9 a 10 12 6 ! Z 23 
Excluded for all foals 3 0 4 4 1 i) 120298) D (22%) 
Sire of all b 4 3 b 4 () 23 (Sbt) 14 (61%) 
9 Sire of some, not all 0 {) : 2 1 { 7 (26%) (27%) 
; BANDS, 2 or more stallions 3 b 2 8 2 6 27 22 
Excluded for all foals 2 2 0 0 ) 2 6 (22%) 4 (18%) 
One 15 sire of al} 1 3 2 4 2 Z 14 (52%) 11 (50%) 
One 15 sire of Sope ( 1 i) 4 0 2 TU E25R) mcts2e) 
FOALS in band with stallion 18 11 23 4) 16 13 121 80 
Without qualifying sire f 3 9 15 i 4 AG (32%) 24 (302) 
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These data suggest that nearly one-third of foals were 
not sired by the harem stallion. The magnitude of this ex- 
clusion rate did not change when sites which may have been 
disturbed by roundup subsequent to conception of foals were 
eliminated from the calculations. The exclusion frequency 
was probably a minimum figure since for some foals no dam was 
identified so blood typing was not at itsS maximum effective- 
ness for detecting incorrectly assigned paternity. 
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Among the 27 harem bands with more than one stallion, 
in only one case did both the dominant and the subordinate 
stallion each appear to sire a foal within the band. Dee i 
possible that in some cases the mares had joined the harem 
subsequent to impregnation by another stallion, although 
Berger (1986) suggests that after recruitment the harem stal- 
lion might force copulation and abort the pregnancy. 


Unfortunately these data do not identify the true sires 
of the excluded foals, whether dominant stallions from other 
bands or stallions from bachelor bands. Not all stallions 
within a populations were captured. Some dominant stallions 
escaped when their bands were being gathered. Without field 
observations from the previous breeding season suggesting 
likely choices, parentage qualification by any of the cap- 
tured stallions could easily lead to spurious conclusions of 
populations dynamics since qualification was taken as proof 
of parentage within bands. 


Band integrity for the months when our populations were 
gathered (September to February) may have been lower than 
during the breeding season and may have been a factor which 
Significantly contributed to the exclusion frequency for 
dominant sires. Field observations during the breeding 
season combined with blood typing would provide the best es- 
timates of breeding by non-dominant stallions or stallions 
from outside the harem band. 


This report corroborates the data of Kaseda, Nozawa and 
Mogi (1982) who combined blood type testing with field ob- 
servations to determine sire-foal relationships among Misaki 
horses which live under conditions similar to wild horses. 
Among 22 foals considered, three were not by the harem stal- 
lion with which the dam associated during the breeding 
season. 


The effectiveness of a locus for parentage determination 
depends on the number of alleles, their frequencies and 
whether genotypes can be directly determined from phenotypes. 
The most effective loci in this study, as in domestic horses, 
were closed blood typing loci (without recessive alleles) 
having five or more alleles with appreciable frequencies, 
namely D, Tf, and Pi. Coat color loci are less effective 
than blood typing for identifying incompatibilities in 
parentage assignment, but could be incorporated as a routine 
check for conclusions from field observations. Among the 
data of this study, five cases of coat color incom- 
patibilities (Trommershausen-Smith, Suzuki & Stormont, 1976) 
were obvious among the exclusions which were detected with 
Dlood typing. Two involved the gray rule, two the chestnut 
rule and one was a palomino without a parent with the coh 
dilution gene. 
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3. Age distribution of qualifying parents of foals 


Using blood typing for verification of parentage, we can 
document ages for the parents of foals, which may be useful 
for predicting dimensions of population increase. These data 
are given in Table VIII. 


Table VIII. Age distribution of qualifying sires and dams 
Trap Site 


F W B Sti €Ai Caz ‘Total (2) 


QUALIFYING SIRE 


N of sires 7 8 7 15 7 = 49 
Average N of foals Lepedie isthe ibe 2.0 © lr4 1.6 
5D Pe Ure Om 0 beAl de br6 0.8 
Range Phe -ZoFl-matssr i-4. 3-2 1-4 
Average age of sire bse tissburthd cid tn ds, 9.1270 l2sd 
SD Si tdzenee pote Hb» 48 a1 
Range Prebipiccom atees Baaze nl, 7°22 ies 


QUALIFYING DAS 


N of daas 15 10 187 miss G2} 17 124 
Average age of daa B.dtebs eee be) Sie] mab 28 Te 
SD awe. wei 0t atid eed 4.0 
Range B16 eialar Salih 257.520, 2226 2-20 


Harem stallions averaged 12.2+5.1 years (range 7-25) and 
dams averaged 7.1+4.0 years (range 2-20), about a generation 
younger than the stallions. These data are strikingly 
Similar to those reported by Perkins et al., 1979 for wild 
horses in the Pryor Mountains of Montana. For harem stal- 
lions Perkins reported a mean age of 12.8+4.0 years (range 
8-20) and for mares a mean of 8.4+3.9 years (range 3-21). 


The likelihood of a mare having a Suckling foal is sig- 
nificantly different among mares grouped by ages (Table IX). 
About one-third of the youngest mares (aged 2-5) and the ol- 
der mares (aged 9-25) qualified to be the dam of a foal, com- 
pared to about two-thirds of those in the 6-8 year age group. 
The difference among the groups is highly significant 
(x2 = 26.78, P<O.0000, 2df). 
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Table IX. Distribution of ages of mares with and without 
foals 


Age of N of With qualifying Without 
mare mares foal 7.) foal (2%) 
2-5 152 ATS CSL) 104 (68.9%) 
6-8 BO Be OOGe ert Soon) 
9-25 oe Zo '€OS.Or:) 49 (66.2%) 
Total SOS 124 (40.7%) 16a 4oOF eo 
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Clearly young mares, including yearlings, can conceive 
and bear a foal, but their ability to raise foals to suckling 
age does not equal that of the 6-8 year old age class. This 
may be due to late onset of maturity in mares at some sites. 
Alternatively, marginal nutritional conditions at some sites 
may not be sufficient to support both growth and reproduction 
in young mares. A young mare with a foal may not conceive 
again for 1-2 years. A mare aged 6-8 is physically mature 
and her own nutritional needs are less, thus she may raise 
foals in successive years. The older mares (9-25) may begin 
to show reproductive decline due to general effects of aging 
or sensitization to blood group factors which can cause foal 
loss due to neonatal isoerythrolysis (NI) 
(Trommershausen-Smith, Stormont & Suzuki, 1975). 


IV. POSSIBLE FUTURE STUDIES 


Although the data of this study appear to provide a con- 
sistent picture of the origins and genetic diversity of wild 
horses of the Great Basin, it cannot be presumed that these 
conclusions would be appropriate for populations of sizes 
less than about 200 animals or for horses outside the Great 
Basin. Thus we suggest that blood typing studies of smaller 
populations and of wild horses at other sites would be 
appropriate. In addition, studies which combine field ob- 
servations during the breeding season with blood typing of 
horses in bands subsequent to foaling would provide essential 
information for interpretation of the paternity exclusions. 
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APPENDIX 1A 
GENE FREQUENCIES 


FLANIGAN 


Nia 


Ls 


BLOOD GROUPS 


A Locus D Locus P Locus 
Aadf: Ose Doc: 0.017 Pa: ees 
Aada: 0.023 Dcgs 0.094 Pb 014 
Aa: 0) DE Oia = O ose 663 
Ab: 0.08S5 Dcgp: O. Ox 
Ac: CHS) Ea | Deowe: @) 
Abc: 0.049 Dcegin: 0.049 Q Locus 
Ae: 0.005 Dcefa: 0. 008 Qabc 0.043 
Abe O Be ge O Qac: ©) 
Ace: 0) Davas O Qb: 0.168 
Abce: 0 DdnOB: Oca ie Ges 0.014 
Aabdg 0) Don7i: O Re O75 
A-: ©. 3538 Ddk: 0.006 
Ddgh: 0.074 
Ddghi9?: 0.014 
C Locus Ddghp: 0) U Locus 
Ca: 0. ba7 Dde: O Ua: 0.445 
Ddeo: 0.06 
Ddeo9B: O 
K Locus Ddek: Ol O1d 
Kas: O Dad: 0.24 
Dadn: O 
Dadé6é9: O 
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